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Introduction 

Due to the prevalence and persistence of cirrus 
cloudiness across the globe, cirrus clouds are believed 
to have an important effect on the climate. Stephens 
1 et al., (1990) among others have shown that the im- 
portant factor determining how cirrus clouds mod- 
ulate the climate is the balance between the albedo 
and emittance effect of the cloud systems. This fac- 
tor was shown to depend in part upon the effective 
sizes of the cirrus cloud particles. Since effective sizes 
of cirrus cloud microphysical distributions are used 
as a basis of parameterizations in climate models, 
it is crucial that the relationships between effective 
sizes and radiative properties be clearly established. 
In this preliminary study, the retrieval of cirrus cloud 
effective sizes are examined using a two dimensional 
radiative transfer model for a cirrus cloud case sam- 
pled during FIRE Cirrus II. The purpose of this pa- 
per is to present preliminary results from the SHSG 
model demonstrating the sensitivity of the bispec- 
tra] relationships of reflected radiances and thus the 
retrieval of effective sizes to phase function and di- 
mensionality. 

Two Dimensional Radiative Transfer Model 

The Spherical Harmonic Spatial Grid (SHSG) 
method of radiative transfer developed by Evans 
(1993) is employed in this study. This method solves 
the two dimensional (2D) radiative transfer equation 
by approximating the radiance field in terms of a 
spherical harmonic expansion for the angular depen- 
dencies and a discrete grid for the spatial properties. 
The discretization of the radiative transfer equation 
in this way, along with the specified boundary con- 
ditions, results in a coupled linear system which is 
relatively sparse and is then solved iteratively using 
the conjugate gradient method.^ The advantage of 
SHSG for this study Is that tbe transfer of radiation 
through arbitrary clouH inhomogeneities can be com- 
puted in an efficient manner. Additionally, since the 
radiance field at all grid points is expanded in terms 
of a spherical harmonic series, the radiances at any 
grid point or angle can be computed by summing over 
the truncated series . The summation of a truncated 


series introduces noise to the solution in the form 
of oscillations. In order to minimize these problems, 
'the 6- AT scaling method is used in order to damp the 
oscillations normally associated with truncating the 
Legendre series of the phase function and the Cesaro 
summation method (Dave and Armstrong, 1974) is 
adapted to two-dimensions to minimize the oscilla- 
tions in the intensities. SHSG also has an indepen- 
dent pixel mode (I.P.) where the horizontal coupling 
of the radiation field is turned off, which simulates 
plane-parallel radiances. The 2D radiances produced 
from SHSG agree with Monte Carlo simulations to 
within 1% in the side scattering portion of the phase 
function, but can differ by at most 8% at the forward 
and backup scattering peaks. At this time the model 
is monochromatic and does not treat molecular scat- 
tering or absorption. 

Estimation of the Cloud Field Optical 
Properties 

For this study a 2D cloud field was obtained from 
a RHI scan by WPL’s Ka-band radar for a cirrus 
cloud sampled on 26 November 1991 during the FIRE 
Cirrus II intensive field observation period (Taneil 
Uttal, personal communication). For the results 
shown here, an image was used which gave cloud re- 
flectivities at 50 meter resolution and represented a 
rather homogeneous looking cirrus cloud character- 
ized by a thick generating cell in the center. The 
cloud reflectivities were converted to ice water con- 
tent using the empirical relation of Sassen (1987), 
1WC = 0.037Z?’ 696 , where 2, is the reflectivity fac-' 
tor of ice (mm 6 ■ m" 3 ) and IWC is the ice water 
content (y * m' 3 ). The estimated ice water contents 
were then used to scale the extinctions computed 
from a modified gamma distribution of equivalent 
area spheres with IWC of 0.0216 g * m“ 3 and an 
effective radius of 80 pm. Table 1 gives the opti- 
cal properties of the cloud simulated in this study. 
The extinctions were computed at the two Landsat 
wavelengths 0.83 and 1.65 pm. The single scatter 
albedos and phase functions were assumed to re- 
main constant across the cloud domain but varied 
with wavelength. The single scatter albedos were 
computed with Mie theory and the phase functions 
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were chosen to have asymmetry parameters consis- 
tent with the geometric optic computations of hexag- 
onal ice crystals by Takano and Liou (1990 and per- 
sonal communication). The asymmetry parameters 
for each wavelength were used as g e /f in the dou- 
ble Henyey-Greenstein function (DHG) to prescribe 
the phase function. The double HG has the form, 
P(coe0) = b gi + (1 - 6) ys, where g u and b are 
constants related by g c /j — b g\ + (\ — b) g 2 . Us- 
ing the preceding assumptions, a cloud field was ob- 
tained which varied in extinction only, which means 
that the effective radius is held constant throughout 
the cloud. 


Table 1: Cirrus cloud optical properties as a func- 
tion of wavelength (A) for the simulations presented 
in this paper. 
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Results 

Normalized nadir reflectances given by R\ = 
(Weilicki et al., 1990) were computed for the cloud 
fields described above using a grid of 1 10x29 which 
represented an area of approximately 11 Arm by 3 
im. The solar zenith angle was 50°. The results are 
shown in the bispectral plot of reflectances at 0.83 
and 1.65 /im as shown in Figure 1. For compari- 
son, a series constant effective radius curves from 15 
fim to 150 /im computed using the phase function 
with 6=1 (DHG1) and the independent pixel ap- 
proximation of SHSG are plotted along with data 
in Fig. 1. As expected, the independent pixel re- 
sults in the top panels of Fig. 1 give relationships 
which agree with the shapes of the effective radius 
curves and the reflectances from DHG1 lie on the 80 
/im curve. However, it is also seen that the other 
forms of the phase function give reflectances that 
give underestimations of effective radii up to 35% for 
$ 0 = 50° and over 80% for 0o = 10°. These under- 
estimations decrease with increasing optical depth 
(optical depth increases from the origin along the 
curves). The sensitivity to the shape of the phase 
function is discussed by Wielicki et ai and is due to 
the magnitude of the phase function at a particular 
scattering geometry. In this case, the bispectral re- 
lationships at f? 0 — 10° are much more sensitive than 
at 0o = 50° because the forms of the phase function 
differed much more greatly in the backscatter por- 
tion (© = 170°) than in the sidescattered portion 


(0 = 130°) of the phase function. The decreased 
sensitivity to the phase function shape as a func- 
tion of optical depth is shown more clearly by Figure 
2 which depicts the retrieved effective radius as a 
function of horizontal position throughout the cloud. 
The solid line gives the ideal retrieval for these cases 
since a distribution of 80 fim was assumed through- 
out the cloud. Referring to the I.P. retrieval lines 
shows that as optical depth increases the difference 
between the curves using DHG1 and DHG2 decreases 
dramatically. This sensitivity to the phase function 
in thinner clouds is due to the single scattering pro- 
cesses which dominate. As the cloud becomes thicker 
multiple scattering becomes more important and the 
individual shapes of the phase function become less 
important to the retrieval of effective radius. 

Figure 1 also shows that unlike the independent 
pixel clouds, the 2D radiances for the upper cluster 
of points do not correspond in slope with the curves 
of effective radius. In fact, the thinner clouds give 
overestimations of effective radius up to a factor of 
2 while the thicker clouds give underestimates of ef- 
fectiv e radius of 10 -25% for Op = 50*. The effect 
of the cloud geometry for this case is to change the 
fundamental relationship between the effective radii 
and the reflectances at the wavelengths shown for 
this preliminary study. This geometry effect can be 
better understood by referring again to Fig. 2. Note 
when comparing the 2D retrieval curves with the I P. 
curves that the differences become greatest for the 
largest optical depth. Referring to the 3rd panel from 
the top and comparing it to the^column optical depth 
in the bottom panel shows that at the front of the 
thickest part of the cloud (i.e., 4.5 km) there are un- 
derestimations of the effective radius while behind 
the cloud there are overestimations. This effect may 
be thought of as an “optical depth” shadowing phe- 
nomena because for this sun angle (50° where the sun 
comes from the left) radiances are enhanced in front 
of the cloud and decreased behind the cloud relative 
to the I.P. calculations. Note that for the higher sun 
angle shown in Fig. 1, this geometric effect is greatly 
reduced. 

Lastly, figure 2 shows a large disagreement be- 
tween the I.P. and curves and the ideal 80 fim line for 
the thinner clouds. This is due to the effects of the 
vertical inhomogeneities. Although this effect is not 
entirely understood at this time, it is clear that for 
thin clouds, as pictured in the top two panels of Fig. 
2, the variation of optical along the vertical has a dis- 
tinct effect on the reflected radiances as compared to 
an homogeneous cloud with an identical integrated 
optical depth. In cirrus which are dominated by thin 
clouds this process may play an extremely important 
role in the interpretation of the radiance fields. 
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Reflection function (1.65 /im) 



Figure 1: Bispectral plots of nadir reflectances for a simulated 2-D cirrus cloud field with variable exinction 
but constant effective radius (r f = 80pm) with the optical properties shown in Table 1. The viewing ange is 
0.0° and <f> c = 180°. The top panels give I.P. results and the bottom panels give 2-D results. The so 


angle is labeled for each column. 
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Conclusions 



Figure 2: Retrieved effective radius as a function 
of horizontal position for the solar geometry of 
e. = 50.0®, <t>* = 0.0® and a viewing angle of 0.0°. 
The top three panels represent 3 different optical 
depths and the bottom panel gives the integrated 
normalized optical depth throughout the cloud. 


These results illustrate that retrieved cirrus effective 
radii contains significant biases not only due to the 
lack of understanding of ice scattering properties but 
also due to the inhomogeneities of the clouds sys- 
tems. For cirrus clouds, which tend to be thinner, the 
sensitivities of the radiance fields to the shape of the 
phase function and to the vertical distribution of the 
optical depth may contribute greatly to biases in the 
retrieval of effective radius. Geometric biases may 
also play a role in very thick cirrus systems which 
sometimes occur in the tropics. Cirrus parameteri- 
zations based upon such retrieved effective radii will 
also produce biases in the radiative properties which 
may be important to the interpretation of climate 
simulations. This study shows that spectral changes 
in the radiance fields of cirrus clouds depend not only 
on the intrinsic microphysical properties of the clouds 
but also upon their spatial structure. 
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